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The generation of entangled photon pairs by parametric down-conversion from solid state CW 
lasers with small coherence time is theoretically and experimentally analyzed. We consider a compact 
and low-cost setup based on a two-crystal scheme with Type-I phase matching. We study the effect 
of the pump coherence time over the entangled state visibility and over the violation of Bell's 
inequality, as a function of the crystals length. The full density matrix is reconstructed by quantum 
tomography. The proposed theoretical model is verified using a purification protocol based on a 
compensation crystal. 



I. INTRODUCTION 

Generation of entanglement is the key ingredient of 
quantum information processing. In optical implementa- 
tion with discrete variables the standard source of entan- 
gled photon pairs is parametric down-conversion in non- 
linear crystals pumped by single-mode laser Recent 
advances in laser diodes technology allow the realization 
of simpler and cheaper apparatuses for the entanglement 
generation 0, H[ , though the quality of the resulting pho- 
ton pairs is degraded by the small coherence time of the 
pump laser. 

In this paper we address theoretically and experimen- 
tally the generation of entanglement using laser diode 
pump as well as its application to visibility and nonlocal- 
ity tests. We focus on the effects of the small coherence 
time and implement a purification protocol based on a 
compensation crystal Q to improve entanglement gener- 
ation. We reconstruct the full density matrix by quan- 
tum tomography and analyze in details the properties of 
the generated state, including purity and visibility, as a 
function of the crystals length and the coherence time 
of the pump. The topics is relevant for applications for 
at least two reasons. On one hand quantifying the de- 
gree of entanglement is of interest in view of large scale 
application. On the other hand a detailed characteriza- 
tion of the generated state allows one to suitably tailor 
entanglement distillation protocols. 

The paper is structured as follows: In Section [IT] we 



describe the experimental apparatus used to generate en- 
tanglement, whereas Section ITTll is devoted to illustrate 
in details the quantum state of the resulting photon pairs 
in the ideal case. The effects of small coherence time are 
analyzed in Section IIVI and the experimental character- 
ization of the generated states is reported in Section [Vl 
Section [VTl is devoted to nonlocality test whereas Section 
IVIII closes the paper with some concluding remarks. 
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II. THE EXPERIMENTAL APPARATUS 

A scheme of the experimental apparatus is shown in 
Fig. [TJ The "state generator" consists of two identi- 
cal BBO crystals, each cut for Type-I down-conversion, 
one half-wave plate (HWP) and one quarter-wave plate 
(QWP) as implemented in ||. The crystals are stacked 
back-to-back, with their axes oriented at 90° with respect 
to each other [l|, Q • The balancing and the phase of the 
entangled states are selected by changing the HWP and 
QWP orientation. 

The crystals are pumped using a 40mW, 405-nm laser 
diode (Newport LQC405-40P), with a spectral line that is 
typically broadened by phonon collisions. The coherence 
time of the pump light t c , which is a fundamental param- 
eter for our experiment, results 544 fs and correspond to 
a spectral width around 0.3 nm. We obtained this im- 
portant information with a standard measurement of the 
first order correlation function. The generated photons 
are analyzed using adjustable QWP, HWP and a polar- 
izer [fjj. Finally light signals are focused into multimode 
fibers which are used to direct the photons to the detec- 
tors. The detectors are home-made single photon count- 
ing modules (SPCM), based on an avalanche photodiode 
operated in Geiger mode with active quenching. For the 
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FIG. 1: Experimental apparatus for generating and analyzing 
entangled states. 



FIG. 2: Geometry for the generation of photon pairs. 



coincidence counting we use a TAC/SCA. 

The nonlinear crystals are properly cut to generate 
photons into a cone of half-opening angle 3.0° with re- 
spect to pump. The first crystal converts horizontally 
polarized pump photons into vertically polarized (V) sig- 
nal and idler photons, while the second crystal converts 
vertically polarized pump photons into horizontally po- 
larized (H) signal and idler photons. This configuration 
introduces a delay time At, depending on the crystal 
length, between the V and the H part of the entangled 
state, as discussed in the following Sections. 



III. THE STATE VECTOR OF THE 
GENERATED ENTANGLED PHOTONS 

The pair of photons generated by SPDC of Type I from 
a single nonlinear crystal, having wave vectors k s and ki , 
are represented by state vectors \k s ) s and \k s )i for the 
signal and idler, respectively. The wavefunction appro- 
priate to the system can be written as a superposition of 
these state vectors @, H, : 



|*) = / d%d% A(u p -^)F(Ak±)f(Ak [ \)\k a ) a \k i )i 

(1) 

where A(u> p — f2°) is the spectral complex amplitude of the 
pump laser, which is a function of the pump frequency 
u) p (kp) = u> s (k s ) + uii(ki), assuming as usual the validity 
of the energy conservation in the generation process, and 
it is centered around the reference frequency 0°. The 
factors F and / are mismatch functions depending on 
the variation of the transverse and longitudinal part of 
the pump wave vector with respect to the reference of 
momentum conservation, and are described in detail in 
the following. 

The function F(Ak±) comes from a spatial integra- 
tion over all the possible processes of photon generation 
within the pump transverse profile in the crystal, taking 
the first order approximation of the nonlinear interaction. 
For a Gaussian pump profile we obtain again a Gaussian 
function, with a width varying as the inverse of the beam 



waist w: 

F{Ak±) = e - w2Ak ^ 4 (2) 

where, referring to Fig. [21 one has 

Ak±_ = k s (uj s ) sin(0 s ) - fci(wj) sin(6>i) 

with the internal generation angles 9 S and 9i for signal 
and idler, respectively. In our case the pump beam waist 
is near 2 mm, therefore we can consider exact transverse 
momentum conservation to a good approximation. In 
fact it is easy to verify that with this beam waist we have 
an angular gaussian width of 0.006° around the reference 
internal angles O s = 0j = 1.8° (derived from external 
angles 9^* = Qf xt = 3.0° using Snell's law), very small 
with respect to the acceptance angle of 0.074° FWHM 
of the optical coupling devices. The conservation of the 
transverse wave vector permits to simplify the geometry 
of the system, by considering in the following a genera- 
tion angle, say 9i, as a function of the other quantities 

LU S , LUi, 9 S . 

The mismatch function /(Afcii) has the same meaning 
of F, but derives from an integration along the crystal 
length Lc, and reads: 



/(Afc,,) = 



sin(Afc||L c /2) 
AkuLc/2 



(3) 



where 



Afc|| = k p (uj p )-k s (u! s )cos(9 s )-k l (u! i )cos[9 i (u! s ,u! il 9 s )} . 

As a matter of fact the pump spectrum width, yet de- 
termining the visibility effects, is very small with respect 
to the spectral width of the down-conversion; this means 
that / is slightly dependent on uj p = lu s + LUi , as can be 
verified numerically. We will not consider such a depen- 
dence by substituting u> p with the reference pump fre- 
quency Clp as the argument of /. This approximation 
turns out to be very good for crystal lengths below a few 
mm, but around 3 mm (our maximum crystal length) 
the conservation of the longitudinal wave vector starts 
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to shrink the down-conversion spectrum. A similar con- 
sideration can be done over the dependence of / over 
the internal angle 8 S ; being the experimental configura- 
tion highly collinear, the optical couplers are practically 
insensible to its variation (within the acceptance cone). 
Therefore we can substitute 9 S with the fixed reference 
angle O s , and the mismatch function becomes: 

f{u v ,w a ,9.) « /(ftp> S) 9 S ) = f(u a ) (4) 

The wavefunction of the photon pair can now be written 
in the simpler form: 

I*) = J dio p doj s dO s A{uj p ~ n°) f(cj s ) 

x \u> s ,6 s ) s \lo p - cj s ,6i(ujp,cj s ,6 s ))i (5) 

In the first approximation we can solve for the integral 
over the internal generation angle 9 S because neither A 
nor / depend on it, but a more refined reasoning put 
forward the fact that the conservation of the transverse 
wave vector introduces a limitation in the effective spec- 
tral width of the mismatch function, hence affecting the 
integration over lo s . This happens because by varying lo s 
around the down-converted reference ftp/2, the idler an- 
gle 9i may go outside from the optical coupler acceptance 
limit, as verified by means of the experimental data dis- 
cussed in Appendix A. This problem does not affect the 
integral over lo v for the smallness of the pump spectral 
width. To take care of this spectral limitation we intro- 
duce a correction factor i?(ft°/2, Acu s ) centered around 
the reference ftp/2 and having the limited spectral width 

Au s (see Appendix A). Defining f(u s ) — f(ui s ) ■ i?, we 
arrive at this wavefunction for the photon pairs: 

I*) = J duj p duj s A(uj p - ftp) f(uj s ) \ui s ) s \uj p - uj s )i (6) 

This expression is used to construct the proper wave- 
function (or the proper state vector) for the entangled 
state generated in our experiment using the pair of ori- 
ented crystals [l|, Q , as described in the previous Section. 
In particular we consider a suitable superposition of the 



single crystal wavefunctions of eq. (J6j) , introducing the 
degree of freedom of polarization on state vectors, be- 
cause the first crystal generates a vertical polarized (VV) 
and the second crystal generates a horizontal polarized 
(HH) photon pairs, respectively. Moreover we have a de- 
lay time between these pairs, due to the different optical 
length of the photon trajectories in the inner of crys- 
tals. This can be represented in the model by assuming 
photon generation in the crystals middle [1(3, HH (in Ap- 
pendix B we show that this is a very good approximation) 
and introducing propagation factors for the internal state 
transport. 

In the Fig. [3] a sketch of the geometry for entangled 
photon generation is shown, limited for clarity to the sig- 
nal trajectories. In the first crystal V photons are gener- 
ated, the state is \V,u> s ) s \V,u> v — u a )i, and the complex 
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FIG. 3: Entangled photon generation and propagation inside 
the crystals. For clarity, only signal photon trajectories (red 
lines) are drawn (idler ones are symmetrically upset). The 
horizontal blue line is the pump ray. (o) and (e) indicates 
ordinary and extraordinary rays, respectively. L c is the length 
of both crystals. 



exponential for the product of the signal and idler prop- 
agation factor (as required by the form of eq. ©) till to 
exiting the crystal is given by: 



P(V) = exp <^ iL 



1 



k?(u p -u) a ) : 



1 



■ k*(w a 



1 



+ k?(uj p - LU S )- 



1 



(7) 



2cos(</>i) ' ' ' o/ 2cos(0i) cos 
where the superscripts (o) and (e) on wave vectors indicates ordinary and extraordinary propagation, and the angles 



2 1 COS{(p 2 )_ 

i rop a 

4>i = 1.807, <f>2 = 1-84 can be found using the laws of wave rays in birefringent crystals |12| and Snell's law, under 
the request of an exit angle of 3°. For the second crystal, in which H photons are generated and the state is 
\H,u) s ) s \H,uj p — bj s )i, the respective propagation factor is: 



P(H) = exp <^ iL c 



k°{uj p ) fc^(wp) 



1 



2cos(0 3 ) 



+ k°(uj p - LU S ) 



1 



2 cos( 



(8) 



where 03 = 1.806, and it has been included the propagation of the pump ray from the generation point of the (VV) 
pair (note that (fis is slightly different from <f>2 due the different refraction index for o and e propagation). The 
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entangled state wavefunction is therefore: 



|tf) = J duo p di^ s A{u Jp ^^l)f{u s )-j={P{H)\H,uj s ) s \H,u Jp ~u Js ) l + P{V)\V,uj s ) s \V,uj p ~u s ) i } (9) 

This expression can be recast in a more useful form in the following way. Let's write the frequencies as io p = tt p + £l p , 
uj s = DP + £1 (with of course £l p = 2fi°), where f2 p and represent the frequency shift with respect to reference for 
the pump and for the down conversion, respectively. Now, in the propagation factors we introduce a first order 
approximation for the wave vectors putting: 

k p (to p ) « k(n° p ) + n p /v p , k s {uj s ) w fc(n°) +n/y, fci(w p - w s ) « fc(Q°) + n p /F - q/f 

where Vp and are the proper group velocities of the pump and of the down converted signal, and these relations must 
be considered both for the ordinary wave and for the extraordinary wave. With these substitutions, and rewriting for 
future convenience the quantum states by factorizing the polarization part from the frequency one, the final form of 
the wavefunction eq. ^ read: 



I*) 



dtt p dQ.A{ti p ) + Q) 



1 

71 



{e i(0H+Tff ^ ) \H) s \H)i\tl) s \n p - Q)i + e «*v-Hvn„) \ V ) s \V)i\n) s \n p - ft),} (10) 



where the phase terms coming from propagation factors are the sum of a constant phase: 



1 " ' ^ pJ ' cos(0 3 ) 
2k°(n°) 4fc e (!!»)U c 
cos(0i) cos(</> 2 ) J 2 



and frequency dependent terms th^ p , Tyf2 p containing 
the total propagation time inside the crystals: 



TH 



TV 



1 

v° 



1 



1 p 



1 



V° cos( 
2 



V° cos( 



V e cos( 



Lc 
2 

Lc 
2 



It is important to note that these delay factors de- 
pend on pump frequency (not on the down converted 
frequency); this fact can be interpreted saying that the 
states (HH) and (VV) exiting the crystals are generated 
from the pump in two different temporal events in the 
past, depending on the different trajectories across the 
crystals. For all these four phase factors, their numerical 
value are determined from the data on refraction indexes 
and group velocities taken from ref. [l3j . and listed in 
the following Table: 





pump 


signal/idler 




(o) 


(e) 


(o) 


(e) 


n 


1.691719 


1.659273 


1.659984 


1.632171 


V 


c/1.77878 


c/1.73901 


c/1.68376 


c/1. 65483 



As a final observation, we note that in writing the fi- 
nal expression for the wavefunction eq. (jlOp . it has been 
discarded the variation of the propagation factors with 
respect to the propagation angles. Due to small angular 
acceptance of the detectors, it is possible to verify that, 
with excellent approximation, this dependence does not 
introduce any relevant effect. 



IV. THE POLARIZATION DENSITY MATRIX 

For the calculation of the density matrix and the com- 
plete characterization of the wavefunction it is impor- 
tant to define at best the statistical properties of the 
CW pump radiation, because our experimental data de- 
pend strongly on its coherence length. In the temporal 
domain, this light is characterized by a (real) constant 
amplitude Aq and a rapidly varying phase with a char- 
acteristic time equal to the coherence time of the pump 
tv.. Therefore we can write: 



du A(u) e % 



A e 



iS(t) 



(11) 



where S(t) is a proper fluctuating phase. The pump am- 
plitude in the temporal domain can be considered as the 
Fourier transform of the complex spectral amplitude over 
a large time interval AT: 



A(u) = — f 
2tt 7 at 



dt Aq e 



(12) 



Our experiment mainly concerns the reconstruction of 
the density matrix of the entangled system on the basis 
composed by the four signal and idler polarization combi- 
nations HH, HV, VH, VV. The relative density operator 
p, from which we derive the reduced density matrix on 
this polarization basis, is obtained from the full density 
operator — \^}(^!\ by tracing over frequencies, e.g. 
by integrating over the frequency state matrix elements: 



p= /'tfc4(ft t > / i (w;-w / | a (w / |*)(*|w / ).K-w / >< (13) 
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corresponding to the fact that we do not perform fre- 
quency measurements. 

The form of the wavefunction in eq. (|10p implies that 
only four elements of the 4x4 reduced density ma- 
trix are different from zero. Using the general relation 
(u)\u> ) = 8(u) — u> ), we straightforwardly obtain for the 
first diagonal element: 



1 



PHH.HH = 



dw|/MI / du p \A(u p )Y 



2 eA ° 2ir 
(14) 

I 



where we put e = f \f(uj)\ 2 and f duj p \A(uj p )\ 2 = 
AqAT/2-k from eq. (fP2|) . With similar calculation 
the other nonzero diagonal element results pvv.vv = 
Phh.hh, as expected by symmetry arguments. 



For the two off diagonal elements one has phh.vv 
Pvv hh ' an d i n particular: 



Phh.vv 



2 / <M/MI 2 



e e 



duip \A(u! p )\ e 
** / du p \A{uo p )\ 2 e- l ^ {TH - Tv) 



(15) 



where we put (f> = (\>h — 4>v ■ With the Wiener-Khinchinc theorem this frequency integral can be recast as a two time 
correlation function over the interval AT, which can be taken very large with respect to the coherence time of the 
pump, and smaller than the detector response time: 



du p \A(uj p )\' 



,(th—tv) 



, AT 



1 



2tt \ATj AT 



rv)) 



A 2 



AT 



-At/t c 



(16) 



r 



where At = \tr — Ty\, and the result is taken from Ref. 
[l|. If Ar»r c we have an incoherent superposition of 
random phases and the average of the complex exponen- 
tials tends to zero, otherwise we have a coherent sum, 
and the integral tends to one. 

Finally, setting the state generator QWP in order to 
have = (see Ref. [l|) and putting for simplicity p = 
e _Ar / Tc , the reduced density matrix is: 





HH 


HV 


VH 


vv 


HH 


( 1 










HV 














VH 














vv 










1 > 



(17) 



which can also be conveniently derived from a sum of two 
distinct density matrix, one of a pure entangled state and 
the other of a statistical mixture: p = p p e + (1 — p)p m 
where p e = |* B )(* e | with |* B ) = -%{\HH) + \VV)) 

and p m = \\HH){HH\ + \\VV){VV\. This is the more 
suitable form used for a comparison between theory and 
experimental data. 



V. EXPERIMENTAL TOMOGRAPHIC 
RECONSTRUCTION OF THE DENSITY 
MATRIX AND CORRELATION VISIBILITY 

In order to fully characterize the generated states at 
the quantum level wc employ quantum tomography of 
their density matrices [l5l ]. The experimental procedure 
goes as follows: upon measuring a set of independent 
two-qubit projectors P M = |'0 At )(V' At | (p = 1, ...,16) cor- 
responding to different combinations of polarizers and 
phase-shifters, the density matrix may be reconstructed 
as Q = iLfiPn^v- where p M = Tv[g P^] are the proba- 
bilities of getting a count when measuring and T p 
the corresponding dual basis, i.e. the set of operators 
satisfying Tr[P M r„] = 8^ [16]. Of course in the ex- 
perimental reconstruction the probabilities p p are sub- 
stituted by their experimental samples i.e. the frequen- 
cies of counts obtained when measuring P M . In order to 
minimize the effects of fluctuations and avoid non physi- 
cal results we use maximum-likelihood reconstruction of 
two-qubit states 0, 0]. At first we write the density 
matrix in the form 



= Ttf, 



(18) 



which automatically guarantees that g is positive and 
Hermitian. The remaining condition of unit trace Tig — 
1 will be taken into account using the method of La- 
grange multipliers. In order to achieve the minimal 
parametrization, we assume that T is a complex lower tri- 
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angular matrix, with real elements on the diagonal. This 
form of T is motivated by the Cholesky decomposition 
known in numerical analysis [l8[ for arbitrary non neg- 
ative Hermitian matrix. For an M-dimensional Hilbcrt 
space, the number of real parameters in the matrix T 
is M + 2M(M - l)/2 = M 2 , which equals the number 
of independent real parameters for a Hermitian matrix. 
This confirms that our parametrization is minimal, up to 
the unit trace condition. 

In numerical calculations, it is convenient to replace 
the likelihood functional by its natural logarithm, which 
of course does not change the location of the maximum. 
Thus the function subjected to numerical maximization 
is given by 



N 

L(f) = In Tr(f f f P^ k ) - ATr(f f T) , (19) 
fe=i 

where A is a Lagrange multiplier accounting for normal- 
ization of q that equals the total number of measurements 
N . This may be easily proved upon writing g in terms of 
its eigenvectors |<^) as g = I] /i y 2 ,|^)(0 M |, with real y^, 
the maximum likelihood condition dLjdy v — reads 



N 



k=l 



yu{(t>u\P^ k \<i>v) 



(20) 



which, after multiplication by y v and summation over u, 
yields X — N. 

The above formulation of the maximization problem 
allows one to apply standard numerical procedures for 
searching the maximum over the M 2 real parameters of 
the matrix T. The examples presented below use the 
downhill simplex or the simulated annealing methods 
[Toj . Results of the reconstruction are reported for crys- 
tals with three different thicknesses, precisely 0.5, 1 e 3 
mm, and in the case of compensation of the delay time 
between generated photons, as discussed later. Moreover 
we present an analysis on the direct measurement of the 
visibility of the entangled state. 



ization angle at 45°. The theoretical prediction is: 
P(£.,45?) = ,(45°| .<£.|p|&>.|45°), 



-p (cos(£ s 



45°)f + -(l-p) (21) 



where £ s is the angle of the signal polarizer in the 
counter-clockwise direction, with the horizontal axis as 
the 0° reference. As it is apparent from this formula, 
when p is near the unity (delay time smaller with re- 
spect the coherence time of the pump) the oscillating 
contribute due to the non-local correlations is dominant. 
On the contrary, with greater delay time (and small p) 
the correlations are washed out and the result is that of a 
statistical mixture which does not depend on the angle. 
In particular, the maximum of P(£ s , 45°) is at 45°, while 
the minimum is at 135°, hence we can write explicitly 
the visibility V of the oscillation as: 



V 



P(45°,45°)-P(135°,45°) 
P(45°,45°) + P(135°,45°) 



(22) 
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0.77 


0.72 


1 mm 


0.50 


0.52 


3 mm 


0.09 


14 



FIG. 5: Tomographic reconstruction of the generated state 
for three different crystals. The measured and calculated vis- 
ibility are shown in table. 




FIG. 4: Entangled state visibility as a function of the polarizer 
angle, for generating crystals of 0.5, 1 e 3 mm thickness. 

Data on correlation visibility are simply obtained by 
removing the HWP and QWP plates of the tomographic 
analyzer (see Fig. [T]) and detecting the signal and idler 
coincidence counts in a time interval, as a function of the 
signal polarizer angle, and having fixed the idler polar- 



fa) (b) 




FIG. 6: Tomographic reconstruction with a delay compensa- 
tion crystal (see text), (a) Crystal angle set for maximum 
compensation, visibility 0.66. (b) Crystal angle at 90° with 
respect to (a), visibility 0.17. 
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In Fig. [4] we show the visibility measurements as a 
function of the signal polarizer angle for the three dif- 
ferent crystal pairs, with the theoretical prediction of 
eq. (f2"2")l indicated by a full line. The comparison be- 
tween the theoretical density matrix elements of eq. (|17p 
and their tomographic reconstruction from experimental 
data is shown in Fig. (5) It is confirmed that the off di- 
agonal elements tend to reduce in magnitude for larger 
crystal thickness; in particular for 3 mm crystals we ob- 
tain the density matrix of a statistical mixture. 

In our model the lack of visibility of the entangled state 
is fully ascribed to the decoherence effect due to the fluc- 
tuating phase difference between H and V parts of the 
SPDC, depending on the delay At. Having a very small 
area of the fiber collimator, we have neglected any deco- 
herence of spatial origin, which introduces a phase varia- 
tion depending on the detector viewing angle. In order to 
verify this statement, we have performed a series of mea- 
surements with the 3 mm crystal, putting windows of 
0.5 mm linear aperture in front of the collimators: if the 
decoherence had a spatial contribution, we would have 
expected an increasing in the state purity. In fact, the 
results of the state reconstruction were the same as the 
original configuration, thus supporting our hypothesis. 

This fact also suggests how to improve the purity of 
the entangled state by a phase retardation on the H po- 
larized part of the pump with respect to the V polarized 
part, to get Sn(t + At) = <5y(i). This can be approxi- 
matively accomplished by inserting, along the pump ray 
and before the state generator, a properly oriented BBO 
crystal with a suitable length. We performed a series 
of measurements using the 1 mm double crystal as state 
generator, and a 3 mm single crystal as pump phase re- 
tarder. By varying the orientation of the axis of this crys- 
tal, we have compensation or enhancement of the effect 
of the time delay between the parts of the generated en- 
tangled state. In particular, the visibility is expected to 
vary from a maximum to a minimum for a 90° change in 
orientation, as confirmed by the tomographic reconstruc- 
tion shown in Fig. [6l Notice that the maximum visibility 
of 0.66 is larger than the corresponding visibility without 
the auxiliary crystal (see Fig. [5]), thus demonstrating a 
partial time delay compensation. 

VI. MEASUREMENTS ON THE VIOLATION 
OF BELL'S INEQUALITY 

We have also performed a series of measurements of the 
S parameter, characterizing the Bell's inequality in the 
CHSH version 20] , for a comparison with the prediction 
of our theoretical model. To obtain reliable data on the 
S parameter we used the same experimental apparatus 
previously employed for correlation measurements. We 
considered as usual the 16 different configuration of the 
polarization angles on the signal an and on the idler (2j. 

The Bell S parameter is theoretically defined as: 

S = E(a, b) - E(a, b') + E(a', b) + E(a', b') (23) 



where the arguments a, a' and b, b' are the selected angles 
for signal polarizer and idler polarizer, respectively. The 
function E is defined as E(a, 0) = P(a, 0) + P(a ± , 0^)~ 
P{a,f3- L )~P{a ± ,f3),wheye a 1 - = a+90° e 1 - = /3+90°. 
The function P is exactly that described in eq. [21] but 
with the idler angle specified by the argument. For any 
realistic local theory one has |iS| < 2, while for quantum 
mechanics | *S* | can be greater than 2, reaching a maximum 
value of 2v2. The following choice for the angles is used: 
a = 0°, b = 9, a' = b + 6 and b' = a' + 6. In this way S 
is a function of the angle 8 alone. In Fig. [7] we show the 
calculated S(9) for three states with different visibility 
V = p = (1,0.7,0.5); by decreasing p, the values of S 
tend to return in the limit of a local theory. 




20 40 60 80 



e n 

FIG. 7: S parameter as a function of 9 for three different 
visibility values. Full line: V = 1; dashed line; V = 0.7; 
dotted line V = 0.5. 
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FIG. 8: Experimental results for S parameter for three values 
of 8, compared with the theoretical S for a visibility of 0.77. 

For a comparison with these results of our model, we 
have measured the S parameter for three different angles 
using as state generator the pair of 0.5 mm crystals, that 
is the case with higher visibility. In Fig. [8] we show the 
theoretical curve of the S parameter for a visibility equal 
to 0.77 (full line), together with two other curves (dashed 
lines) indicating the extremal of the experimental errors, 



8 



relative to the limited number of count during data ac- 
quisition. The three measurements of S for the angles of 
16°, 24°, 40° are indicated with error bars. In particular 
we get S{16°) = 2.38 ± 0.03, S{2A°) = 2.417 ± 0.025 and 
5'(40 o ) = 0.80 ± 0.05. From these data we can conclude 
that in the case of 24° the Bell's inequality is violated for 
more than 17 standard deviations. 



VII. CONCLUSIONS 

We have analyzed, both theoretically and experimen- 
tally, the generation of polarization-entangled photon 
pairs by parametric down-conversion from solid state 
CW lasers with small coherence time. In particular, we 
have analyzed in some details a compact and low-cost 
setup based on a two-crystal scheme with Type-I phase 
matching. The effect of pump coherence time on the 
entanglement and the nonlocality has been studied as a 
function of the crystals length. The full density matrix 
has been reconstructed by quantum tomography and the 
proposed theoretical model is verified using a purification 
protocol based on a compensation crystal. We conclude 
that laser diodes technology is of interest in view of large 
scale application and that its that the characterization of 
the generated state allows one to suitably tailor entan- 
glement distillation protocols. 
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APPENDIX A: MEASUREMENT OF THE 
COHERENCE LENGTH 

In this Appendix we experimentally verify that the 
spectrum of the down-converted signal is reduced when 
coincidence photon counts are performed, as a conse- 
quence of the trasverse momentum conservation. 

If we observe only a single photon of the generated 
pair, the part of the spectrum incident on the coupling 
device is described in practice by the mismatch function 
f{uj s ) defined in Eqs. © and (HI). But if we observe both 
photons and measure the simultaneous counts between 
signal and idler, we will detect a spectrum with a smaller 
width, and therefore we have a greater coherence length 
of the radiation. This because if we have a very wide 
spectrum for the signal at a fixed angle of observation, 
the idler photons, correlated with the signal photons also 
by transverse momentum conservation, will be dispersed 
over an angle that can be wider with respect to the ac- 
ceptance angle of the coupling device. Hence the pair 
of coupling devices work as a filters limiting the spectral 



window for observation. To the purpose of a determina- 
tion of this effective spectral width, we present here some 
measurements using interference methods. In particular 
we performed two series of measurements, the first rela- 
tive to the direct counts on a single detector to find the 
width associated with f(u s ), the second relative to the 
coincidence counts on the two detectors to determine the 
width of the corrected mismatch function f(u> s ) used in 
eq. ©. 




Interferometer 



FIG. 9: Sketch of the experimental apparatus for the measure 
of the coherence length. 

In Fig. [9] we show the experimental scheme (based on 
a single BBO crystal) employed for these types of mea- 
surements. An interferometer equal to that described in 
Q is placed among the signal ray. This interferometer is 
easy to align and is a very stable device. In both types 
of counting measurements we expect to see interference 
fringes as a function of the delay time introduced by the 
interferometer between two optical paths, and within the 
radiation coherence time. In particular we would deter- 
mine a greater coherence length in the case of coincidence 
counts with respect to the case of signal single counts. 

The theoretical description of the interferometric ex- 
periment is as follows. In the case of a single photon ob- 
servation, and with a crystal generating horizontal pho- 
tons, the density matrix for the signal before the inter- 
ferometer can be built with the wavefunction of eq. (|10|) 
of Section IIIII by neglecting vertical polarization states, 
and tracing over the idler frequency: 

pi = J dn\f(n° + n)\ 2 \H,n) s s {H,n\ (ai) 

where we do not consider the immaterial propagation 
factor and use the original mismatch function of eq. 

The density matrix after the interferometer follows by 
considering: (a) a polarization rotation of 45° due to the 
HWP plate placed before the first calcite crystal; (b) the 
delay time r introduced by the interferometer between 
the H and V parts; (c) the projection of these states 
over the axis of the final polarizer oriented at 45°, placed 
before the coupling device. The final density matrix is 
easy obtained as: 
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pi = J dn\f(n°+n)\ 2 ^ |i + e *° T | 2 \45°,n) s s {45°,n\ 

(A2) 

The probability to observe a count on the detector is then 
proportional to: 



Pi(r)i 



dn s (45°,n |pi|45°,n 



dfi\f(n° + n)\ 2 - |i + e ' 



(A3) 



The width of the interference pattern representing count 
numbers as a function of the delay t, is given by a fac- 
tor similar to a Fourier transform of the down-converted 
power spectrum; hence this width scales as the inverse of 
the spectral power width of the function /. 

In the case of signal and idler coincidence counting, the 
state vector is again derived from eq. (jTUJ) , by taking only 
the H part and discarding the propagation factor. After 
the passage in the interferometer, it is straightforward to 
see that the coincidence probability is the same as for 
the single count probability by replacing /(f2° + fi) with 
the modified mismatch function /(fi° + f2). Hence in this 
case the width of the interference curve is governed by 
the modified spectral power width Aw s of /. 



(a) 



(b) 
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FIG. 10: Interference patterns: (left) single signal counts; 
(right) signal and idler coincidence counts. 

In Fig. [10] we show on the left the interference pattern 
obtained with signal single counts, using the BBO crystal 
of 3 mm length. The width of the curve is near 30 fs, 
corresponding to a down converted spectrum of about 64 
nm. On the right we show the pattern in the case of 
signal and idler coincidence counts: the coherence time 
is enlarged to 70 fs, corresponding to a spectral width of 
27 nm. In both cases the coherence length is well below 
that of the pump light. These data are used to determine 
the appropriate correction factor R(Auj s ) — /(w s )//(oj s ) 
in the definition of the wavefunction eq. ([B]) . 



die. But in fact these states can be generated in any point 
in the inner of the crystals, therefore the propagation fac- 
tors P(H) and P(V) must be position dependent. Let's 
indicate with z\ and z 2 the longitudinal coordinates of 
the internal generation points for the first crystal and for 
the second crystal, respectively. Referring to the Fig. [3j 
we now have the following two equations for the propa- 
gation times th and Ty\ 



th{zi,z 2 ) 



Lc - zi 

V° 
1 p 



L 



c-zi 



V° costs' 



(Bl) 



y°cos(0i) V e cos(0 2 ) 



(B2) 



Generally speaking, the state visibility p would depend 
on the delay time At(zi,z 2 ) = th(zi,z 2 ) — tv{zi,z 2 ). 
Because we do not have any information about the effec- 
tive position in which a particular photon pair is gener- 
ated, we consider an average over the possible positions, 
by integrating with a flat distribution probability: 



Lc' 



dz\ dz 2 e 



-At(zi,z 2 )/t 



(B3) 



In Fig. [TJJ we show a comparison between the visibil- 
ity for state generation in the crystals middle, and that 
obtained from the above formula. It is clear that there 
is some difference only for very small visibilities obtained 
with very long crystals. 




APPENDIX B: COMPLETE CALCULATION OF 
THE DELAY TIMES IN STATE GENERATION 

In deriving the delay time between (HH) and (VV) 
photons, we assumed state generation in the crystals mid- 



FIG. 11: Entangled state visibility as a function of the crystal 
length for the different assumptions on the state generation 
position. Full line: in the crystals middle; Dotted line: in the 
whole crystal length. 
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